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PACS 05 . 70 . Ln - Nonequilibrium and irreversible thermodynamics 

Abstract - We apply the model of minimally nonlinear irreversible heat engines developed by 
Izumida and Okuda [EPL 97, 10004 (2012)] to refrigerators. The model assumes extended Onsager 
relations including a new nonlinear term accounting for dissipation effects. The bounds for the 
optimized regime under an appropriate figure of merit are analyzed and successfully compared with 
those obtained previously for low dissipation Carnot refrigerators in the finite-time framework. 
Explicit derivation of the Onsager coefficients is also presented. 



Nowadays, the optimization of thermal heat devices is 
receiving a special attention because of its straightforward 
relation with the depletion of energy resources and the 
concerns of sustainable development. A number of dif- 
ferent performance regimes based on different figures of 
merit have been considered with special emphasis in the 
analysis of possible universal and unified features. Among 
them, the efficiency at maximum power (EMP) ?7maxP for 
heat engines working along cycles or in steady-state is 
largely the issue most studied independently of the ther- 
mal device nature (macroscopic, stochastic or quantum) 
and/or the model characteristics [TI-[T4]. Even theoretical 
results have been faced with an experimental realization 
for micrometre-sized stochastic heat engines performing a 
Stirling cycle [T5] . 

For most of the Carnot-like heat engine models ana- 
lyzed in the finite-time thermodynamics (FTT) frame- 
work [316] the EMP regime allows for valuable and simple 
expressions of the optimized efficiency, which under en- 
doreversible assumptions {i.e., all considered irreversibili- 
ties coming from the couplings between the working sys- 
tem and the external heat reservoirs through linear heat 
transfer laws) recover the paradigmatic Curzon-Ahlborn 
value [E] ?7maxP = 1 - = r]GA using t = Tc/Th, where 
Tc and Th denote the temperature of the cold and hot heat 
reservoirs, respectively. A significant conceptual advance 
was reported by Esposito et al. [17j by considering a low 



dissipation Carnot heat engine model. In this model the 
entropy generation in each heat exchange process is as- 
sumed to be inversely proportional to the time duration 
of the process, the reversible regime is approached in the 
limit of infinite times and the maximum power regime al- 
lows us to recover the Curzon-Ahlborn value 77CA when 
symmetric dissipation is considered, but without the re- 
quirement of assuming any specific heat transfer law nor 
the linear-response regime, since the derivation is indepen- 
dent of the external temperature values. Considering ex- 
tremely asymmetric dissipation limits, these authors also 
derived that EMP is bounded from the lower and upper 
sides as 



-77- S ?7maxP S -X : 

2 2 - 77c 



(1) 



where rjQ = 1 — t is the Carnot efficiency. The Curzon- 
Ahlborn value 77CA [i.e.., the symmetric dissipation limit) 
is located between these bounds. Additional studies on 
various heat engine models P"5^Bn] confirmed and gener- 
alized above results. 

Work by de Tomas et al. [21] extended the low dissipa- 
tion Carnot heat engine model to refrigerators and besides 
proposed a unified figure of merit (^-criterion described 
below) focusing the attention in the common character- 
istic of every heat energy converter (the working systems 
and total cycle time) instead of any specific coupling to 
external heat sources which can vary according to a par- 
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ticular arrangement. For refrigerators, we denote the heat 
absorbed by the working system from the cold heat reser- 
voir by Qc, the heat evacuated from the working system 
to the hot heat reservoir by Qh and the work input to the 
working system hy W = — Qc- For heat engines, we 
may inversely choose natural directions of heat transfers 
and work. That is, we denote the heat absorbed by the 
working system from the hot heat reservoir by Qh, the 
heat evacuated from the working system to the cold heat 
reservoir by Qc and the work output by = Qh — Qc. 
Then we can define this ^-criterion as the product of the 
converter efficiency z times the heat Qin exchanged be- 
tween the working system and the heat reservoir, divided 
by the time duration of cycle icycio: 



X 



zQi, 



cycle 



(2) 



where Qin = Qc, z = e = Qc/W for refrigerators and 
Qin = Qh, z = 1] = W/Qh for heat engines. It becomes 
power as X = W^/icycio when applied to heat engines and 
also allows us to obtain the optimized coefficient of per- 



formance (COP) 



-maxx 



under symmetric dissipation con- 



ditions when applied to refrigerators: 



1 



;CA- 



(3) 



This result could be considered as the genuine counterpart 
of the Curzon-Ahlborn efficiency for refrigerators. It was 
first obtained in FTT for Carnot-like refrigerators by Yan 
and Chen |22j taking as target function eQc, where Qc is 
the cooling power of the refrigerator, later and indepen- 
dently by Velasco et al. |23II24| using a maximum per-unit- 
time COP and by AUahverdyan et al. [25] in the classical 
limit of a quantum model with two n-level systems inter- 
acting via a pulsed external field. Very recent results by 
Wang et al. [5B] generalized the previous results for refrig- 
erators HU and they obtained the following bounds of the 
COP considering extremely asymmetric dissipation limits: 



^ S^maxx — 



V9 + 8ec~'S 



(4) 



where ec = ^jiX — r) is the Carnot COP. These bounds 
for refrigerators play the same role that eq. ([T]) for heat 
engines. 

A shortcoming of the all above results is the model de- 
pendence. Thus additional research work has been de- 
voted P7H2^ to analyze its validity and generality using 
the well founded formalism of linear irreversible thermo- 
dynamics (LIT) for both cyclic and steady-state models 
including explicit calculations of the Onsager coefficients 
using molecular kinetic theory [30fl32] . Beyond the linear 
regime a further improvement was reported by Izumida 
and Okuda |33j by proposing a model of minimally non- 
linear irreversible heat engines described by extended On- 
sager relations with nonlinear terms accounting for power 
dissipation. The validity of the theory was checked by 




J2 



Fig. 1: Schematic illustration of the minimally nonlinear irre- 
versible refrigerator described by eqs. (|12[l and (jlSp . 

comparing the results for EMP and its bounds with eq. ([!]). 
Now, the main goal of the present paper is to extend 
the application of this nonlinear irreversible theory |33) 
to refrigerators and analyze their COP under maximum 
X-condition. To get this we first present the model and 
analyze the results when the tight-coupling condition is 
met. Then, we show that the low dissipation Carnot re- 
frigerator model [55] is a particular case and we derive the 
corresponding Onsager coefhcients. 

Refrigerators are generally classified into two types, that 
is, steady-state refrigerators and cyclic ones. Our theory 
can be applied to both cases. Since the refrigerators come 
back to the original state after one-cycle, the entropy pro- 
duction rate a of the refrigerators agrees with the entropy 
increase rate of the heat reservoirs as 



W Qc 



Tc 



W 



I 
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(5) 



where the dot denotes the quantity per unit time for 
steady-state refrigerators and the quantity divided by the 
time duration of cycle tcycic for cyclic refrigerators. From 
the decomposition of & into the sum of the product of ther- 
modynamic flux Ji and its conjugate thermodynamic force 
Xi as a = J2i=i 2 ^i^ij '^^^ define Ji = x, Xi = F/Ty, 
and J2 = Qc, X2 = i/Th — f/Tc for steady-state refrig- 
erators. Then W = Fx where F and x are the time- 
independent external generalized force and its conjugate 
variable, respectively. Likewise we can also define 



Jl = 



I 



^cyclc 
J2 = Qc, X2 = 



Xi = 

1 



w 
I 



(6) 
(7) 



for cyclic refrigerators, where Xi and X2 play the role 
of driver and driven forces, respectively. Anyway, we as- 
sume that the refrigerators are described by the following 
extended Onsager relations |55] : 



Jl — LiiXi 
J2 — L21X1 



L12X2, 
L22X2 



IcJl', 



(8) 
(9) 
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where L^j's are the Onsager coefficients with reciprocity 
Li2 = L21. The nonlinear term —j^jf means the power 
dissipation into the cold heat reservoir and 7c > denotes 
its strength. We also disregard other possible nonlinear 
terms by assuming that they are too small when compared 
with the power dissipation term. 

The Onsager relations are usually written by using only 
the linear terms with respect to Xi and X2 and then the 
non- negativity of the entropy production rate a = JiXi + 
J2X2 leads to restriction to the Onsager coefficients [301 

MM 



Lii > 0, L22 > 0, L11L22 — L12L21 > 0. 



(10) 



Although our model includes the power dissipation as the 
nonlinear term, we also assume eq. (jlOp holds. This idea 
that the dissipation is included in the scope of the linear 
irreversible thermodynamics is also found in . 

The power input P is given hy P = W = JiXiT^ > 0. 
Then the heat flux into the hot heat reservoir (jh = P + 
Qc = J3 is given as 



J3=P + Qc^JlXiTh + J2. 



(11) 



By solving eq. (|S]) with respect to Xi and substituting it 
into eqs. ([9]) and (fTTj) . we can rewrite J2 and J3 by using 
Ji instead of Xi as (see fig. [I]): 

J2 = ^Ji+i22(l-9')^2-7cJi', (12) 
^11 

J3 = J^^Ji+L22il-q^)X2+lhJi^ (13) 

^11 Jc 



where q 



is the usual coupling strength paramc- 



ter (|g| < 1 from eq. (fTO)) ) and 7h > denotes the strength 
of the power dissipation into the hot heat reservoir as 



7h 



ill 



■7c 



(14) 



Eqs. and allow a clearer description of the re- 
frigerator instead of eqs. ^ and ^ by considering Ji as 
the control parameter instead of Xi at constant X2. So, 
each term in eqs. (fT2l) and ([T3|) has explicit physical mean- 
ings (see also |33]): the two first terms mean the reversible 
heat transport between the working system and the heat 
reservoirs; the second ones account for coupling effects 
between the heat reservoirs; and the third ones account 
for the power dissipation (~ J^) into the heat reservoirs, 
which inevitably occurs in a finite-time motion (Ji 7^ 0). 

Using eqs. and ([13]) instead of eqs. dH) and 
the power input P = J3 — J2 and the entropy generation 
a 



^ — are written as 



P 



L21 



1 Ji 



^Ji' > 0, 
^11 



(15) 



and 



In the following we restrict our analysis to the case \q\ = 
1 called the tight-coupling condition. Below, we will show 
that this condition holds for the low symmetric dissipation 
Carnot refrigerator by calculating its Onsager coefficients 
iy 's explicitly. Generalization of the result for jgl 7^ 1 is 
straightforward although cumbersome. 

For the refrigerator, the converter efficiency z in eq. ^ 
is the COP £ defined as 



P 



•II 
P 



ai(i-l)Ji + f-j2 

11 \ T ' Lii J- 



(17) 



where we used eqs. (|12p and We also introduce the X" 
criterion as a product of the COP e and the cooling power 
J2, identifying the cooling power Qin/icycie in cq. ^ as 
Qc — J2- 



J2- yfe-^l-^cJl^ 

p ^(i-l)J, + ^J,- 



(18) 



ill 



where we used eqs. (fT2|) and (fTSl) . Though eq. ^ us- 
ing icycio appears to be valid only for cyclic heat energy 
converters, eq. ([TS]) using J2 is valid even for steady-state 
refrigerators. We consider the optimization of this crite- 
rion with respect to Ji. From eq. (fT2l). we see that Ji 
should be located between and L2i/(7ciii), in order to 
get a positive cooling power as it must be in a refrigerator. 
dx/dJi = reduces to the cubic equation: 



J 3 , 3L2iLii7c — 2L2iScTh , 2 

27cLiiec7h 

9 r 2 t3 
^^21 T , ^21 



-Ji 



0. 



(19) 



Then the sole physically acceptable solution [d'^x/dJi^ < 
0] is given by: 



^1 



2L21 



37ciii + 7ciii JSec (1 + ^) + 9 



(20) 



Substituting eq. (|20p into eq. (HZ]), we obtain general 
expression of the COP under maximum x-condition for 
the model of minimally nonlinear irreversible refrigerators 
with the tight-coupling condition: 



9 + 



See ( 



1 



7h 

7c 



2(1 + ^ 



(21) 



Eq. ((2T)) is the main result of this paper. It is a monotonic 
decreasing function of 7h/7c and we can obtain the lower 
and the upper bounds by considering asymmetric dissipa- 
tion limits 7h/7c — > 00 and 7h/7c ~> 0, respectively, as 



7h 



> 0. 



(16) 



^ ^maxx — 



V9 + Sec - 3 



(22) 
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in full agreement with results reported in |26] and ex- 
pressed in above eq. 

Now we prove that the above model of minimally nonlin- 
ear irreversible refrigerators with the tight-coupling con- 
dition includes the low dissipation Carnot refrigerator [55] 
by explicitly calculating the Onsagcr coefficients. This low 
dissipation Carnot refrigerator is an extension of the qua- 
sistatic Carnot refrigerator by assuming that heat trans- 
fer accompanying finite-time operation in each isothermal 
process is inversely proportional to the duration of the 
process tf. 



In particular, under symmetric conditions (Sh 
further optimization by a gives a 



Sc) the 

maxx = \/l-T+ 1, frOm 



5h = ThA5 



(23) 
(24) 



where we denote by AS* > the quasistatic entropy 
change of the working system during each isothermal pro- 
cess, and by Si/ti the corresponding entropy production 
with a constant strength > [l7l[2Tl[26]. Since it 
works as a cyclic refrigerator and according to eqs. ^ and 
([7|). we can define its thermodynamic fluxes and forces as 
Ji = 1/teycie = !/((« + l)ih), J2 = Qc, ^1 = W/T^ and 
X2 = 1/Th — 1/Tc, respectively, with a = tc/^h- By using 
these definitions and eqs. (|23| and (|24| . we can calculate 
the Onsager coefficients Lij's and the strength of the dis- 
sipation 7i's as 



Lii L12 
L21 L22 




Y 



7h = rhi;h(a + 1), 

rcSc(a + l) 
7c = : 



(25) 

(26) 
(27) 



respectively, where the common factor Y in cq. (|25p is 
given as 



y = (a + 1) ThSh 



TcSr 



(28) 



as expected before, 
L21 and the tight- 
= 1 hold. Thus our 



Prom eq. ((25)) it is obvious that, 
both the Onsager reciprocity L12 

coupling condition lol = 

model eqs. (fT2|) and (fT3)) under the tight-coupling condi- 
tion |g| = 1 includes the low dissipation Carnot refrigera- 
tor in eqs. (f23| and ([24]) as a special case. 

Finally wc note that we obtained eq. (j2ip maximizing x 
only by one-parameter J\ , whereas in the original low dis- 
sipation model the optimization uses the two-parameters 

and [2T1[26]. Thus the parameter space for the max- 
imization is different. In the formalism presented here, 
a further optimization of x = x(q!) by a is still possible 
after the optimization by J\ and it corresponds to the two- 
parameter maximization done in |21[|26j . But, that fur- 
ther optimization does not change the bounds in eq. (j4|. 



which we can straightforwardly obtain the symmetric op- 
timization COP in eq. (3) pT| . 

In summary, we introduced the extended Onsager re- 
lations as the model of minimally nonlinear irreversible 
refrigerators under the tight-coupling condition. We de- 
rived the general expression for the coefficient of perfor- 
mance (COP) under maximum ^-condition. Calculating 
the Onsager coefficients and the strength of the power dis- 
sipation explicitly, we proved that our model includes the 
low dissipation Carnot refrigerator model as a particular 
case. Then, the validity of eqs. ([3]) and ^ derived in the 
low dissipation framework, is confirmed by the minimally 
nonlinear irreversible formalism presented here. 
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